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We propose to induce QCD by fermions in the adjoint representation of the gauge 
group SU (N c ) on the lattice. We consider various types of lattice fermions: chiral, 
Kogut-Susskind and Wilson ones. Using the mean field method we show that a first 
order large- iV phase transition occurs with decreasing fermion mass. We conclude, 
therefore, that adjoint fermions induce QCD. We draw the same conclusion for the 
adjoint scalar or fermion models at large number of flavors Nf when they induce 
a single-plaquette lattice gauge theory. We find an exact strong coupling solution 
for the adjoint fermion model and show it is quite similar to that for the Kazakov- 
Migdal model with the quadratic potential. We discuss the possibility for the adjoint 
fermion model to be solvable as N c — > 00 in the weak coupling region where the 
Wilson loops obey normal area law. 
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1 Introduction 



This paper is inspired by a recent work by Kazakov and Migdal []J who proposed to 
induce QCD by a scalar field on the lattice. This scalar field is taken in the adjoint 
representation of the gauge group SU(N C ) while the gauge field is attached in the usual 
way to make the model gauge invariant except no kinetic term for the gauge field. This 
construction makes the model solvable in the large- N c limit. The latter circumstance 
differs the Kazakov-Migdal model from the previous attempts to induce Wilson lattice 
gauge theory by fundamental-representation scalars or fermions H at the number of 
flavors (or different species) Np ~ N c . 

A price for the solubility is an extra local symmetry which has to be spontaneously 
broken |4| under the way to the continuum limit. As we proposed in Ref. [||], if the large- 
N phase transition occurs for the Kazakov-Migdal model before the one separating the 
confinement and perturbative Higgs phases, one gets normal confinement (area law) in 
the intermediate region bounded by the two phase transitions. The problem is, however, 
that this scenario in not realized for the original model of Ref. |TJ . 

The mean field analysis || indicates that no first order phase transition occurs, with 
decreasing mass parameter, before the gaussian model becomes unstable which is com- 
patible with the exact solution of Refs. |6|, [7| obtained solving the master field equation 
and reproduced by means of loop equations ||. The instability is due to an unlimited 
condensation of scalar particles for an action which is unbounded from below and were 
lead to the Higgs phenomenon if a stabilizing self-interaction would be added. 

The latter fact prompts the following way out. Let us consider fermions rather than 
bosons. On the one hand, the induced lattice action obtained by the large mass expansion 
looks for fermions very similar to that for scalars. On the other hand, since fermions 
never condense, the gaussian fermionic action in an external field is stable so one can 
arbitrarily decrease the mass parameter. Moreover, if one assumes that the gauge field 
is weakly fluctuating, the perturbative expansion for fermions looks like that for QCD 
and we expect, therefore, a phase transition which separates the perturbative region from 
the strong coupling phase with unbroken Zn symmetry and associated local confinement. 
This phase transition is to be identified with the large- N phase transition discussed above. 

In the present paper we propose to induce QCD a la Kazakov-Migdal by fermions in 
the adjoint representation of the gauge group on the lattice. We consider various types 
of lattice fermions: chiral, Kogut-Susskind and Wilson ones. We discuss first the ad- 
joint scalar and fermionic models at large number of flavors Nf and show they induce a 
single-plaquette lattice gauge theory quite similarly to Refs. [0, [| for the fundamental 
representation. Using the mean field method we show that a first order large-iV phase 
transition occurs for chiral and Kogut-Susskind fermions at Nf = 1. We conclude, there- 
fore, these two models induce QCD. The similar conclusion presumably holds for Wilson 
fermions while the results for Nf = 1 are less certain. We find an exact strong coupling 
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solution of loop equations for the adjoint fermion model and show that for chiral fermions 
it is quite similar to the solution [[jj for the Kazakov-Migdal model with the quadratic 
potential while for Wilson fermions the exact result coincides with the leading order of 
the large mass expansion. We point out that, as usual, the strong coupling solution is 
not sensitive to the first order large- N phase transition which occurs inside the region 
where the large mass expansion is convergent. We discuss the possibility for the proposed 
models to be solvable as N —>■ oo in the weak coupling region where the Wilson loop 
averages obey normal area law. 



2 Inducing single-plaquette action at large Nf 

The scenario of breaking the Zn discussed in Section 1 can be completely realized at 
large number of flavors Nf when the induced gauge theory is in our case the single- 
plaquette lattice gauge theory with the adjoint action quite similarly to Refs. p], Q for 
the fundamental representation of SU(N C ). The large- N phase transition occurs for the 
former action according to the pattern we proposed long time ago ||. 



2.1 Scalars 

For simplicity let us begin with the case of scalars. The extension of the Kazakov-Migdal 
model IJ to the case of Nf flavors reads 

where Nf fields (i = l,...,Nf) take values in the adjoint representation of the 

gauge group SU(N C ) and the link variable U^(x) belongs to the gauge group. For the 
purposes of the present paper, we restrict ourselves by the quadratic potential. The 
original Kazakov-Migdal model corresponds to Nf = 1. 

Calculating Nf gaussian integrals, one rewrites the model ( [2.1Q as the lattice gauge 
theory 

Z = f HdU^x) e - s - Mx)] , (2.2) 

with the induced action given by the following large mass expansion 

N f \trU(T)\* 

i>ind[U\- — ^ 2J(r) • 

* r /(r)r?V 
For Nf = 1 one recovers the result of Ref . PJ . 

We apply now the idea by Bander [[J of how to obtain the single-plaquette lattice 
action from the sum over loops on the r.h.s. of Eq. ( |2.3| ). Let Nf amd m\ simultaneously 
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tend to oo so that ttIq ~ (Nf)&. Then the only single-plaquette term survives while those 

_i 

associated with longer loops are suppressed at least as (Nf) 2 . Therefore we have shown 
in this limit that 

S md = -^- s J2\trU(dp)\ 2 ( 2 - 4 ) 

where dp is the boundary of the plaquette p. In other words the induced theory is just 
the single-plaquette lattice gauge theory with the adjoint action. 

As is known |§ the lattice gauge theory with the action (|2.4j ) undergoes the (first 
order) large- N phase transition at Nf/(Amfy ~ 2. This is exactly the phase transition 
advocated for the Kazakov-Migdal model by Kogan, Semenoff and Weiss [§j] and needed 
to induce QCD according to the scenario of Ref. 0. Below the phase transition, the 
adjoint action ( |2.4| ) is reduced || at N c = 00 to a fundamental single-plaquette action 
and the Wilson loop averages obey area law. 

It is crucial for this consideration that the critical value 

1 

m *~(l^) 4>jD (2 ' 5) 

as Nf — > 00 so that the gaussian model ( |2.1|) is stable. A mechanism to escape the 
instability while having the large- N phase transition is pretty simple: the contribution of 
long loops in the effective action (|2.3|) which were lead to condensation is now suppressed 
by construction. This pattern looks like the idea to avoid instability by means of fermions 
since in the fermionic case the long loops do not condense due to alternating signs in the 
sum (see Eq. (|2.9|) below). 

On the contrary there is no such phase transition for m\ > D at Nf = 1. This 
statement relies on the mean field study of Ref. || and is compatible with the exact solu- 
tions 0) HI- Therefore, we expect the large- N phase transition to occur for the model ( |2.1| ) 
only for Nf larger than some value Nf*. This value is estimated in Section 3.5 by the use 
of the mean field method. 



2.2 Fermions 

Similarly to how Hamber || used large-iV/ limit to induce the Wilson single-plaquette 
lattice gauge theory by fermions in the fundamental representation of SU(N C ), we can 
extend the results of the previous section to the fermionic case. 

Let us consider the adjoint fermion model (AFM) which is defined by the partition 
function 

z AFM = /n^wnn^w^/W e " SFM,!/ ' ^) 

x,p. X f=\ 
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Here Sf[$, ^f, U] is the standard action for fermions on the lattice: 

N f 

S F [*,V,U] = ]T]TiV c tr (mVf(x)* f (x) 

f=l X 

-J2^ f (x)p-U,(x)^ f (x + i,)Ul(x) +* f (x + riP+Ul(x)* f (x)U,(x)}) (2.7) 

where ^ f(x) is the Grassmann anticommuting field while 

Pi = r ± 1, (2.8) 

stand for the projectors. The case r = corresponds to chiral fermions while r = 1 is 
associated with Wilson fermions. As is well known, the chiral fermions describe 2 D iVj 
flavors in the naive continuum limit while Wilson fermions are associated with Nf flavors. 

The path integral over \1/ can be calculated by the large mass expansion similarly to 
the scalar case. One gets for the induced action 

where Sp stands for the trace over the spinor indices of the path-ordered product of the 
projectors P^ (plus or minus depends on the orientation of the link /) along the loop V. 

The contribution of leading order in 1/m comes from a plaquette which gives 

^- 2f ^ 2r2 - r4) £W)|-. (2.io) 



m P 



If m ~ Nj , only this term survives in the large- Nf limit while the others are suppressed 
_i 

at least as Nf 2 similarly to the scalar case. 

The estimate of the critical values of m which are associated with the large- iV phase 

i i 
transition gives m (2iVj)4 for chiral fermions and m* ps (47V/) 4 for Wilson fermions. 

As we shall see in the next section, the mean field analysis indicates that the phase 

transition certainly exists in the former case at Nf — 1 as well. 



3 Mean field analysis 

While we have shown in the previous section that both adjoint scalar and adjoint fermion 
models induce the single-plaquette Wilson lattice gauge theory in the limit of large Nf, 
there is no hope for this limit to be solvable even at iV c = oo (this were correspond to 
an exact solution of Wilson lattice gauge theory at iV c = oo). We investigate in this 
section AFM with Nf = 1 which, as is discussed in Section 5, could be as solvable as the 
Kazakov-Migdal model. 
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3.1 Mean field with fermions 



We apply to AFM defined by the partition function ( |2.6| ) at Nj = 1 the (variational) 



mean field approach of Ref. ||. Since we substitute the mean field solely for the gauge 
field and exactly integrate over fermions, we expect the mean field method is applicable 
for our case despite the fermionic system. 

Substituting the mean field value [U^(x)]ij = j]8y for all links of the lattice except the 
given one, we get the following self-consistency condition 

2^-^) = ! (3.1) 



where 

/ Ux d$?(x)dV(x) e -SH*.*.»?U tr (*(0)P"*(0 + fi) + tf(0 + //)P+# (0 



'A 



and 



(3.2) 



- V 2 Y / [^(x)P~^(x + ^ + *(x + M )P+^x)]) . (3.3) 

Eq. ( p.l|) holds providing rj > 1/2, which is prescribed by the requirement for the cor- 
responding one-matrix model to possess a weak coupling solution. To find the location 
of the phase transition one uses either geometric criterion which is based on terminating 
the 77 7^ solution with increasing m or the thermodynamic criterion which says that the 
phase transition occurs when free energies of the 77 = and rj 7^ phases coincide. 



3.2 Chiral fermions 

For chiral fermions the gaussian integral on the r.h.s. of Eq. ( |3.2|) can easily be calculated 



by the standard lattice technique (see Ref. HlPp to give 



(3.4) 

Here 1q{ol) is the modified Bessel function. 

The r.h.s. of Eq. (|3.4j) multiplied by if monotonically increases with decreasing m 



approaching the maximal value 2 D /D at m — 0. Moreover, it can be approximated with 

' m 2 
•2r) 4 



2 

great accuracy by the first term of the expansion in x = l/(f^T + D) which can be seen 



from the series expansion 

rfb A = 16x - 8x 2 + 32x 3 - 42x 4 + 168x 5 - 320x 6 + 1280x 7 - + 

1397319^12 , 1007910^13 , <n( 



-x w + 124047a; 11 - i^a; 12 + 1397319x" + 0(x 14 ). (3.5) 
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Thus, we substitute 

^2- D {^ + DW (3-6) 

and there is no need in numerical calculations. 

The critical value m* can be easily determined from the above mentioned geometric 
criterion which says 

drj 



dm 



dm 

oo or — — 
drj 



(3.7) 



at the critical point m*. Differentiating Eq. with the r.h.s. substituted by Eq. 
one gets the solution 



2 

m* « 6.5, and 77* ~ - (3i 



at £> = 4. 



Some comments about the obtained critical values of the large- N phase transitions in 
the case of chiral fermions are in order. 

• In contrast to the scalar case analyzed in Ref. 0, the integral on the r.h.s. of 
Eq. (O) never diverges for < m 2 < oo which is a consequence of the stability of 



the fermionic system discussed above. Since i] 2 bA can now reach values larger than 
0.5, Eq. ( |3.1| ) always possesses a solution. 

The coefficient 2 D on the r.h.s. of Eq. ( |3.4|) is nothing but the number of flavors which 

is associated with the naive continuum limit of chiral lattice fermions. Therefore, 

one could relate the appearance of the large-iV phase transition with an effective 

number of flavors similarly to Section 2.2. This point of view explains why, as is 

shown below, the large- N phase transition still certainly occurs for Kogut-Susskind 

d ' 
fermions describing 2 2 flavors and presumably for Wilson lattice fermions which 

correspond to 1 flavor in the naive continuum limit. 



3.3 Kogut— Susskind fermions 

As is well known, the number of continuum flavors described by chiral fermions can 
be reduced by 2"? passing to Kogut-Susskind fermions. To this aim one performs the 
following transformation of the spinor $(x) 

y(x) = ( 11 y^ l2 y\..( lD y° x (x) (3.9) 

where the coordinates of the lattice vector x^ are measured in the lattice units. In the 
new variables the fermionic lattice action ( |2.7| ) reads 

Sf\x,X,U) =J2 N ° tT ( m x(x)x(x) 

X 

D 

~ E ^(x)[x(x)U^x)x(x + v)Ul(x) - X (x + fi)Ul(x)x(x)U,(x)]) (3.10) 
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where 

771(2:) = 1; -q^x) = (3.11) 

Since the action ( |3.10| ) is diagonal w.r.t. spinor indices, one can keep only one component 
of xi x ) so that it can be considered as a scalar Grassmann variable. The normal spinors 
are reproduced in the continuum limit by 2 D x's taken at neighbor sites of the lattice. 

D 

Such a procedure manifestly reduces the number of flavors in the continuum by 2 2 . 
The mean field analysis of AFM with Kogut-Susskind fermions is quite similar to that 

D 

of the previous section. One should only replace 2 D in Eqs. (|3.4j ) end ( |3.6|) by 2 2 so that 
the self-consistency condition (|3.1| ) at D = 4 reads 

2 (--l) =^ + 1. (3.12) 

This equation determines the critical values associated with the large- N phase transition 
to be 

m* « 0.7, and 77* w -. (3.13) 

Therefore, we have shown that the large-iV phase transition occurs for AFM with 
Kogut-Susskind fermions. The critical values (|3.13j ) are smaller than those given 
by Eq. (|3.8|) and are near the limit of applicability. 



3.4 Wilson fermions 



For AFM with Wilson fermions b& is given by Eq. (|3.2|) . The standard calculation [10] of 
the gaussian fermionic integral yields 



22 " r* dp p 1 



in ^-E M cosp M 



1--^ ^ ~\ ^ I. (3.14) 

V 2r] (^-Z.cosp,) +E M sin 2 p M/ 



The integral on the r.h.s. of Eq. (|3.14 ) possesses a number of interesting properties. 
Its large mass expansion 

fbA ^(m\ m (m\o((^r) (3,5) 

begins with the term ~ m -4 rather than m~ 2 as for chiral fermions or bosons. The point 
is that backtracking paths in the sum over paths associated with 6^ are now forbidden 
since the product of projectors 

P+P- = r 2 -l (3.16) 

vanishes for r — 1. For this reason our mean field coincides for Wilson fermions with the 
one which would be obtained directly from the induced action ( |2.9|) where the backtrack- 
ing paths never contribute due to unitarity of U . A price for this nice property is the 
complicated integral on the r.h.s. of Eq. ( p,14| ). 
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Since b A varies between at m = oo and 1 at m = 0, one expects Eqs. (|3.1|) and (|3.7|) 
to possess a solution for the critical value m*. We performed an approximate numerical 
analysis and found m* ~ 77* « 0.5 which looks like the case of Kogut-Susskind fermions. 
However, now an uncertainty of the results is bigger so that our conclusion is that, while 
the large-iV phase transition presumably occurs for Wilson fermions, more accurate in- 
vestigation maybe by other methods is needed. However, the large- iV phase transition 
certainly occurs for AFM with Wilson fermions for Nf > 2. 

3.5 Nf scalars 

The mean field method can be applied to estimate the critical value Nf* below which 
the large- N phase transition discussed in Section 2.1 disappears. The corresponding self- 
consistency condition is given by Eq. (|3.1| ) where b A reads 

/Vf r°° i^q. 
b A =-± dae'^ a l$-\a) L(a) (3.17) 
2r] z Jo 

and 1(a) being the modified Bessel function. 

The condition rj > 1/2, which is prescribed by the requirement for the corresponding 
one-matrix model to possess a weak coupling solution, yields rfbA > 1/2 or Nf* > 30 - 
a pretty big value — since the maximal value of the integral on the r.h.s. of Eq. Q3.17Q is 
0.06 at D = 4 . 



4 Loop averages at strong coupling 

AFM can be exactly solved at N c = 00 in the strong coupling limit (i.e. before the 
large- N phase transition). We obtain the result solving loop equations quite similarly 
to Ref. || where the Kazakov-Migdal model with the quadratic potential was solved by 
this method. The solution for AFM with chiral fermions obtained below is quite similar 
to that by Gross |7| for the Kazakov-Migdal model with the quadratic potential and 



coincides with the old results by Klugberg- Stern et al. [fn]] for the infinite-coupling limit 
of lattice N c = 00 QCD with chiral fermions. We show this fact is not a coincidence and 
is related to very peculiar properties of Wilson loop averages in AFM model as well as 
in the Kazakov-Migdal model. For Wilson fermions the exact solution of loop equations 
coincides for the same reason with the leading order of the large mass expansion. It will 
be explicitly demonstrated this strong coupling solutions are not singular, as usual, at the 
point of the first order large- iV phase transition which is determined by another criterion. 
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4.1 Loop equations in AFM 



All the observables in AFM can be expressed via the adjoint Wilson loop 

W A (C) = ^(\trU(C)\ 2 -l)^ , (4.1) 

where the average is understood with the same measure as in Eq. (|2.6| ), and the fermionic 
path-dependent amplitude 

Gij(C xy ) = (i-tr (y t (x)U(C xy )VMU\C xy )) ) (4.2) 

where i and j are spinor indices. 

The amplitude Gij(C xy ) obeys two sets of loop equations which are analogous to the 
loop equations in QCD with quarks (in the fundamental representation of the gauge 
group). It is simple to obtain an equation resulting from the Dirac equation for ^>(x) 
which reads 

D 

mG(C xy ) - £ [P+G(C {x+ ^ )x C xy ) + p-G(C (x ^ )x C xy )} = 5 xy W A (C xy ) (4.3) 

where the matrix multiplication over spinor indices is implied. Here the path C( x+fl ) x C xx 
is obtained by attaching the link (x, /i) to the path C xx at the end point x as is depicted 
in Fig. 1. The projectors are given by Eq. (|2.8|) for chiral and Wilson fermions or by 
Eq. (|3.11|) for Kogut-Susskind fermions. Due to the presence of the delta-function, Wa 
on the r.h.s. of Eq. (fO|) always enters for closed loops, i.e. it is an object of the type ( fll 



Similarly to Ref. [§] where an analog of Eq. (fO|) was solved for the Kazakov-Migdal 
model with the quadratic potential at strong coupling, Eq. (|4.3j) admits the ansatz 




G t ,(C xy ) = G L \ J] Pt\ (4-4) 



where plus or minus corresponds to the direction of the link / which belongs to the contour 
C xy and L is the algebraic length |12| {i.e. the one after contracting the backtrackings) of 



C xy . The spin factor is needed to cancel the projectors in Eq. ( fl.3| ). 



For such an ansatz, Eq. ( ^4.3| ) for AFM is reduced to the following recurrent relations 
on Gl 

mG L - G L _i - (2D - l)oG L+1 = for L > 1 , 

mG - 2DaGi = 1 (4.5) 

where a = r 2 — 1. These equations are quite similar to the corresponding equations for 
the Kazakov-Migdal model with the quadratic potential which are given || by the same 
formulas with a = 1 and m replaced by 2ml. The drastic simplification of the loop 
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equation (|4.3|) is due to the fact that Wa(C) vanishes in the strong coupling region at 
N = oo except for C with vanishing minimal area A m i n (C) (i.e. contractable to a point 
owing to unitarity of U's): 

W A (C) = So,A min (o) + O (^) . (4.6) 
Eq. (|4.5| ) can be solved introducing the generating function 

oo 

G(X) = E G L \ L . (4.7) 

L=0 

The solution at o ^ (the solution for a = i.e. for Wilson fermions is described in 
the next section) is expressed via Go which was up to now arbitrary. It can be deter- 
mined imposing || the single-pole analytic structure of G(A) as a function of the spectral 
parameter A. This requirement unambiguously determines G(X) to be 

= 2(2 C -l) ffG „ 

{^m? + 4(1 - 2D)a - m) A + 2{2D - I) a 

and 

C = ^ (4.9) 

m(D - 1) + Dwm 2 + 4(1 - 2D)a 



The single-pole prescription is justified in the next section. 



4.2 Properties of strong coupling solution 

Let us consider the strong coupling solution of AFM, obtained in the previous section, 



from the viewpoint of the path representation of the amplitude (12). At N c = oo one 

gets 

GiAC X y) = E J^W A (C xy T yx ) ( I] If) (4.10) 

where the path C xy Y yx is closed by construction. Eq. ( |4.10|) is nothing but the standard 
large mass representation JTO] of the amplitude in lattice gauge theories with fermions 



while the Wilson loop is taken in the same representation as fermions (adjoint in our 
case). The fact that the only one- loop average enters Eq. ( |4.10| ) is due to N c = oo 



In the strong coupling region where Wa{C) is given by the simple formula (|4.6p , the 
r.h.s. of Eq. ( |4. 1U[ ) drastically simplifies. The path T yx must coincide with C xy passing 



in the opposite direction modulo backtrackings^ since Wa equals 1 in this case and 
otherwise. Therefore, the problem of calculating the amplitude G by the sum over paths 

2 One should not confuse these backtrackings of T with possible backtrackings of C discussed above. 
The latter ones can be always contracted due to unitarity. 
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is reduced to a problem of counting backtrackings of one-dimensional trees embedded in 
a D-dimensional hypercubic lattice. The above formula ( |4.8| ) represents the solution of 
this problem. Notice that since the spin structure on the r.h.s. of Eq. ( [4.10| ) factorizes, 



the coefficient functions Gl are quite similar for chiral fermions and scalars. 

We would like to point out that exactly the same problem emerges in lattice QCD at 
infinite coupling where quarks are in the fundamental representation but the Wilson loops 
for non-contractable paths vanish since the coupling is infinite. Its solution obtained for 
chiral fermions by Klugberg-Stern et al. fill] using a different method exactly coincides 
with Eq. (pi). 



The result of Ref. [J 1 was reproduced by Mkrtchyan and one of the authors [13] by 
means of loop equations. While the loop equation resulting from the Dirac equation for 
fermions looks like Eq. ( |4.3| ) with Wa on the r.h.s. substituted by the fundamental Wilson 
loop, the loop equations resulting from the (lattice) Maxwell equation are simplified since 
fermions are now in the fundamental representation. All the definitions in this case are 
the same as above with the path ordered phase factor in the adjoint representation substi- 
tuted by the one in the fundamental representation, say the definition of the fundamental 
amplitude G reads 

Fin \ — /x,F(„\TT(rt \tf,F 



G^(C xy ) = (x)U(C xy )^ („)). (4.11) 

The equation resulting from the variation of U^(x) at the link (x, fi) G C xy emanating 
from the point x reads[| 

G F (C xy ) = G F (C x(x+ ^h ll G F (C xy ) - G F (0)^G F (C {x+ ^ y ) (4.12) 



for chiral fermions. For the ansatz ( [4.4[) one gets 



G L = GqGl-i - oGxG L for L > 1. (4.13) 

While the corresponding loop equation for AFM does not coincide, generally speaking, 
with Eq. ( |4.12| ), for the ansatz (fOj) it should reduce, as is discussed above, to Eq. ( |4.13| ). 
At L — 1 this equation gives one more relation between Go and Gi which together with 
the L = equation (|4.5| ) result in a quadratic equation for Go having an unambiguous 
solution ( |4.9| ) expandable in 1/m . Therefore, we have shown that the above strong 
coupling solution for AFM with chiral fermions (or for the Kazakov-Migdal model with the 
quadratic potential) is unique. In particular, we have justified the single-pole procedure 
of Ref. §. 

An analogous exact strong coupling solution of loop equations for AFM with Wilson 
fermions can be obtained as well. Since a = for r = 1, Eqs. (|4.5| ) are further simplified 
so that the solution reads 

Gl = —t (4-14) 



3 The bilinear in G terms were missing in the original paper by Weingarten 
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which coincides with the leading order of the large mass expansion. The reason behind 
this it as follows. As is discussed in Section 3.4, the problem of backtrackings is missing 
for Wilson fermions due to the property ( |3.16 ) of projectors P^. Therefore, the leading 
term of the large mass expansion gives the exact result for AFM with Wilson fermions. 

Some comments about the properties of the exact strong coupling solutions ( |4.8| ) and 
( [4.14j ) of AFM with chiral and Wilson fermions are now in order: 

• The solutions (fL8|), fl£I3p are not singular for m 2 > in contrast to the scalar 
case. Therefore, the strong coupling solution is fermion case is stable everywhere. 
However, AFM undergoes the first order large-iV phase transition after which the 
strong coupling solution is not applicable. As usual for first order phase transition, 
this critical point is quite regular from the viewpoint of the large mass expansion and 
is determined by another method (like the mean field approach of Section 3). After 
the large- N phase transition occurs, Eq. ( |4.6| ) that corresponds to local confinement 
is no longer valid and should be replaced by normal area law. The ansatz ( O ) is no 
longer applicable in this region and the solution of loop equations would be quite 
different. 

• While for D = 1 the strong coupling solution ( |4.8| ) coincides with the free theory 
amplitudes, it looks simpler for D > 1 than the corresponding free theory results 
given by formulas which are similar to those of Section 3. 

• All the statements of Ref. concerning the breaking of chiral symmetry for QCD 
with chiral fermions at infinite coupling are applicable to our case. In particular, 
no extra phase transition associated with spontaneous breaking of chiral symmetry 
occurs with decreasing mass and < > which is determined by Eq. ( [4.9[ ) is 
non- vanishing at m — 0. 

• Loop equations turned out to be perfect to solve the backtracking problem arising 
in the large mass expansion and unambiguously determine the solution without any 
assumptions. 



5 Conclusions and discussion 

The main conclusion we draw from the results obtained in the present paper is that a sim- 
ple replacement of scalars by fermions in the Kazakov-Migdal model cures the problem of 
instability and provides the large- N phase transition when the fermion mass is decreased. 
While in the strong coupling phase (before the phase transition) one gets local confine- 
ment which is characterized by Eq. fl4.6p , normal area law holds in the weak coupling 
phase (after the phase transition). This phase resembles very much the standard lattice 
gauge theory. Moreover, at large- Nf a single plaquette lattice action arises. Hence, we 
conclude that AFM induces QCD. 

The obtained exact strong coupling solution of AFM does not know anything about 
the large- N phase transitions which is associated with freezing of the gauge field U^{x) 
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near some mean field value r\ while 77 = in the strong coupling region. This is usual for 
first order phase transitions. The weak coupling solution of AFM should be quite different 
from the strong coupling one. 

After a nonvanishing value of 77 emerges, the arguments of Ref. |I[ about the continuum 
limit caused by long loops become applicable. We expect therefore that AFM will induce 
continuum QCD as well. One more argument in favor of this is the mentioned above 
similarity between the weak coupling phases of AFM and of the Wilson lattice gauge 
theory. 

A very interesting question is whether the solubility of the Kazakov-Migdal model 
by means of the master field equation || would be pursued for AFM. This were give a 
way to solve AFM in the weak coupling region which would be an alternative to loop 
equations. As is mentioned already, we do not think this would be the case for any Nf. 
We see however some hope to derive an extension of the Itzykson-Zuber-Metha integral, 
which played a crucial role in the approach of Refs. ^ , to the case of Kogut-Susskind 
like fermions: 

where ip and \ are Grassmann variables which transform under the adjoint representation 
of SU(N C ). While it is not clear of how to reduce ip and x to a diagonal form, the integral 
over the Haar measure on the r.h.s. of Eq. (|5.1| ) could be tractable by other methods. 
This problem deserves future investigations. 



Note added 

When this paper was being prepared for publication, there appeared two more papers 



I~5| , |l6l| on the Kazakov-Migdal model. While Ref. |l5j deals mostly with the quadratic po 



tential, a strong coupling solution for a general model with fundamental Wislon fermions 
is obtained by Migdal . 
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Fig. 1 The graphic representation for Gij(C xy ) (a) and Gij{C^ x+ ^ x C xy ) (b) entering 
Eq. ( f4.3|) . The empty and filled circles represent *$?i(x) (or *$?i(x + //)) and *&j(y), 
respectively. The oriented solid lines represent the path-ordered products U{C xy ) 
and U{C { x+ ^) x C xy ). The color indices are contracted according to the arrows while 
the spinor ones are associated with circles. 
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